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Decoherence suppression via environment preparation 
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To protect a quantum system from decoherence due to interaction with its environment, we 
investigate the existence of initial states of the environment allowing for decoherence-free evolution 
of the system. For models in which a two-state system interacts with a dynamical environment, we 
prove that such states exist if and only if the interaction and self-evolution Hamiltonians share an 
eigenstate. If decoherence by state preparation is not possible, we show that initial states minimizing 
decoherence result from a delicate compromise between the environment and interaction dynamics. 

PACS numbers: 03.65Yz, 03.67-a, 03.67Pp 



Introduction Decoherence provides an elegant frame- 
work to explain why an open quantum system coupled 
to its environment will exhibit a set of preferred states, 
usually ruling out a coherent superposition of arbitrary 
states. In Zurek showed that even if the superposition 
principle treats all quantum states equally, the interac- 
tion between the quantum system and its environment 
would select a restricted number of "pointer" states (ein- 
selection) and destroy the phase coherence of superpo- 
sitions of those pointer states (decoherence). This phe- 
nomenon presents a formidable challenge for such appli- 
cations as quantum computation. In this simplest model 
of decoherence, it was readily realized that initial states 
of the environment exist that allow for decoherence-free 
unitary evolution of the quantum system. These pecu- 
liar states were usually neglected on the basis that "in 
realistic cases, such highly ordered initial states [. . . ] are 
unlikely to be relevant" |2j or because the environment 
self-evolution would preclude such a unitary evolution 

In this paper, we investigate the conditions under 
which such initial states of the environment do exist in a 
framework where the quantum system interacts with its 
environment and the environment also evolves by itself. 
The results obtained underline the crucial role of the en- 
vironment's self-evolution. The ability to identify and 
prepare such special initial states could be used in order 
to store quantum states. Indeed, even if the environment 
dynamics cannot be controlled, it might be possible to 
prepare it in a specific initial state. However, our results 
restrict what can be expected from such a technique. 

More precisely, we obtain a mathematical condition for 
the existence of an initial state allowing decoherence-free 
evolution in the presence of an interaction Hamiltonian 
and a self-evolution of the environment, stated in terms of 
the structure of the two Hamiltonians. Next, we analyze 
in detail a particular model, which is an extension of the 
model introduced by Zurek in Finally, we assess the 
impact of imperfect state preparation and discuss how to 
choose an initial state that minimizes decoherence when 



state preparation cannot avoid it altogether. 

The model We consider a quantum system S and 
its environment £ with dynamics described by a Hamil- 
tonian of the form 



H = S<S>H+1®H £ . 



(1) 



The Hermitian operator S acts in the Hilbert space of 
<S, which we take to be two-dimensional. The quantum 
system is thus taken to be a quantum bit (qubit) Q. the 
Hermitian operators H and Hg act in the Hilbert space 
of the environment. 

The total Hamiltonian (l} induces pure dephasing and 
is typical of a coupling between the system and the en- 
vironment that commutes with the self-evolution of the 
system (jj. 

Without loss of generality, we can assume S — a z . Its 
eigenstates, of eigenvalues ±1, are written |0) and |1), 
respectively. 

Suppose that the global system is in a product state 
at t = 0: 

|*(o)> = |V(o)> ® \i) 

where \tp(0)) — a |0) + b\l) is an arbitrary normalized 
pure qubit state. At a later time t, the state evolves to 

|*(t)) = o|0)®|e (t)) + 6|l><8)|ei(t)> 

where 



i—\e k (t))=H k \e k (t)), 



0,1, 



(2) 



and where we have defined Hq = Hg + H and Hi = 
Hs — H with initial condition |eo(0)) = |ei(0)) = \I). 

Since \^{t)) is no longer a product state in general, the 
reduced density matrix p(t) of the quantum system S no 
longer describes a pure state: the system has decohered. 
To quantify this, the off-diagonal elements of p(t) in the 
basis {|0) , |1)} are reduced by a factor 

r(t) = (e (f)|ei(*)>. 
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Essentially, r(t) quantifies how distinguishable the states 
\so(t)) and \si(t)} are. If they differ only by a phase, i.e. 
\r(t)\ = 1, they are indistinguishable and the system S 
remains in (or has returned to) a pure state. If r(t) = 0, 
they can be distinguished with certainty, decoherence is 
complete, and the state of the quantum system reduces 
to a statistical mixture \a\ 2 |0) (0| + |6| 2 |1) (1|. In general, 
the system has only undergone partial decoherence and 
a straightforward calculation shows that its purity is 

P(t)=Trp 2 (t) = l-2\a\ 2 \b\ 2 (l-\r(t)\ 2 ) 

where |f*(i)| 2 is the Loschmidt echo @|. 

Thus, a qubit prepared in an arbitrary coherent super- 
position of pointer states will remain pure if and only if 
the environment can be prepared in a state |J) for which 

\r(t)\ = 1 for alU (3) 

i.e., for which \so(t)) and \£i(t)} are the same state up to 
a time-dependent phase at all time. 

Condition for coherent evolution Given that the state 
evolves according to ^ , what can be said of the Hamil- 
tonians H and Hs in order that for some initial state |J), 
\so(t)) and \si{t)) satisfy the the coherence criterion ([3])? 
The factor r(t) is given by 

r(t) = {I\e tHot e- tHlt \I} . 

Clearly, if \I) is an eigenstate of both Hq and Hi, then 
r(t) is a phase and ([3]) is satisfied, so the existence of 
a common eigenstate of these Hamiltonians (or, equiva- 
lent^, of H and Hs) is a sufficient condition. It is also 
a necessary condition, although the state \I) itself need 
not be a common eigenstate. To understand this, write 
H in terms of its spectral decomposition, with eigenval- 
ues and associated projection operators {nj -*} (so 

that Ho = Ej A j° )n j 0) )> and similarly for Hi. Then 

r W = ^e^ 0, -^)* (/ |nfn«|/}. (4) 

In general, r(t) contains terms of different frequencies. 
However, in order for it to remain of unit magnitude for 
all t, only one frequency can appear; all others must be 
associated with vanishing coefficients. The easiest way 
for this to occur is if \I) is a common eigenstate of Hq 
and H\ , in which case the double sum collapses to a sin- 
gle term. A more general possibility is if \I) is a linear 
combination of common eigenstates of Ho and Hi corre- 
sponding to the same energy difference. In terms of the 
original Hamiltonians H and Hs, \I) must be a linear 
combination of eigenstates of Hs and each of these must 
also be a degenerate eigenstate of H. Stated otherwise, 
we must be able to find a basis in which H — Cm ffi H' 
and Hg = Dm®H' £ , where Cm is an M -dimensional con- 
stant matrix (proportional to the identity) and Dm is an 



M-dimensional diagonal matrix; |/) can be any vector in 
the first M dimensions. 

Thus, preparing an initial state of the environment per- 
fectly maintaining the coherence of the qubit at all times 
can only be done if both the interaction Hamiltonian and 
the environment self-evolution Hamiltonian exhibit a spe- 
cific structure, namely they share an eigenstate. 

Can we characterize the rarity of such pairs among 
all pairs of Hamiltonians? From a mathematical point 
of view, it can be shown that pairs of Hermitian ma- 
trices sharing an eigenstate are a closed set with empty 
interior. Intuitively, this implies that pairs of Hermitian 
matrices with a common eigenstate are rare. This conclu- 
sion is also supported by the observation that even if both 
Hamiltonians share a common eigenstate, most perturba- 
tions will destroy this property. To see this, suppose that 
the Hamiltonians are written in a basis where H is diago- 
nal and suppose that only the first eigenstate is common 
to both Hamiltonians. Then the first row and column of 
Hg are zero except for the (1,1) element. In order for a 
perturbation of Hs to preserve the common eigenstate, it 
(the perturbation) must have the same structure, so the 
real dimension of perturbations preserving this structure 
is less than the real dimension of all possible perturba- 
tions ((AT - 1) 2 + 1 vs A 2 ). 

Notice that if the environment self-evolution Hamilto- 
nian Hg is zero, any eigenstate of the interaction Hamil- 
tonian is a suitable initial state. Thus, it is the dynamics 
of the environment that usually prevents the existence 
of such an initial state. In the following section, we will 
make this situation explicit by computing the loss of co- 
herence induced by adding dynamics to an otherwise- 
static environment. 

Environment evolution as a perturbation In this sec- 
tion, we consider a solvable model due to Zurek [l( to 
which we add a new term 1 eg) Hs to provide dynam- 
ics to an otherwise-static environment. We consider the 
case in which the environment is initially prepared in an 
eigenstate of H and show that the self-evolution of the 
environment will destroy the coherence of the system. In 
this model, the environment consists of n spin-1/2 par- 
ticles. The dimension of the environment Hilbert space 
is N — 2™. Every particle of the environment interacts 
with the system through a a z a z interaction. The eigen- 
states of the Pauli matrix a 2 , acting on the fc th spin are 
{|0) fe , |l) fc }- The strength of the interaction between the 
system and the k th spin of the environment is measured 
by a coupling constant gk € R. Thus, the Hamiltonian is 

n 

H = a z s ®Y,9^l (5) 

fe=i 

The eigenstates of H = J2k=i 9*,®% are the states in 
which the fc th spin of the environment is in an eigen- 
state of cr^, i.e., the states \x) — (££)j. =1 \xk)k where all 
.Xfe G {0, 1}. A convenient way to order these states is to 
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consider that a; is a number between and N — 1 of which 
the binary representation is x = x\X2 ■ ■ ■ x n . Thus, 

Va;e{0,l}" H\x)=w x \x) 

where uj x — J2k=i(~^) Xk 9 k - Throughout this section, we 
will assume that the coupling constants {gk} are chosen 
so that the eigenvalues u> x are distinct. 




Figure 1: Decoherence model with dynamical environment 

So far, the spins of the environment do not interact 
with each other; thus, preparing the environment in any 
state \x), being an eigenstate of H, will give rise to co- 
herent evolution of the system. We will now add a self- 
evolution of the environment of the form 



H E = A 



E 

z,J/6{0,l}™ 



x) (y\ 



where A e R is a perturbation parameter. This 
Hamiltonian is proportional to the projector on |$) = 

T^EzefO,!}" \ x ) = ®fe=i l+fc) (+fcl> where \+k) = 



V2 



|0*> 



7fl lfc > : 



H £ = XN |$) ($| 



Clearly, the interaction and the environment Hamiltoni- 
ans do not share a common eigenstate because all eigen- 
states of H have a small overlap with |<f>). 

Suppose that the environment is prepared in the state 



10} 



'10), 



(6) 



fc=i 



which is an eigenstate of H with eigenvalue ujq. 

Let us compute the decoherence factor r(t) for this 
initial state of the environment. Standard perturbation 
theory for A <C min x ye i iy n \oj x — ui y \ shows that 



Ht)\ 



sin 4 (E^EoA 



The time average of this quantity is given by 

1 



|r(t)| 2 = 1-6A 2 



y . 

^ (CJ - U x ) 2 



Therefore, the time average is reduced by a factor propor- 
tional to J2 x ^o (u a -u )' 2 > wm ch characterizes the density 
of energy levels near the unperturbed energy level. 



Imperfect preparation of the environment Suppose 
now that the condition for coherent evolution in the orig- 
inal model (HJ is satisfied, so that there is at least one 
state for which the time evolution preserves the coherence 
of the system S. We have assumed implicitly our ability 
to prepare perfectly the initial state of the environment. 
Notice that this strong assumption is nonetheless weaker 
than requiring control over the environment dynamics 
at all times. The interest of preparing the environment 
has already been studied from both a theoretical and 
an experimental Q point of view. However, preparing 
the environment in a given initial state is a difficult task 
which might only be achieved partially. Let us give a 
simple example in the case of Zurek's model ©. When 
trying to prepare the state iJH), the goal is to prepare 
all spins in the state ak |0} fc + (3k with ctk = 1 and 
(3k = 0. Suppose that we are only able to ensure that 
\(3k\ 2 < £ <§C 1 for all k, i.e., all spins of the environment 
are prepared with a small error. In that case, the average 
value |r(i)| 2 is bounded by 



|r(t)| 2 > ((l- £ ) 2 + £ 2 ) 



2\" £<1 



1 - 2ne 



which is attained if \(3k 



e for all k. Thus, for small 



independent errors on each of the n spin of the environ- 
ment, the coherence loss is proportional to n. This par- 
ticular example indicates that even if an initial state al- 
lowing for decoherence-free evolution exists, the difficulty 
to prepare it will grow with the size of the environment, 
as one would expect intuitively. 

Imperfect control of the environment We now address 
a situation which is in a sense opposite to the one just 
considered. Rather than having perfect control over the 
dynamics (so that a common eigenstate of H and He, can 
be made to exist) but an imperfect ability to prepare the 
initial state, suppose that we can prepare perfectly any 
state we wish but that no such common eigenstate ex- 
ists. Given that the coherence of S cannot be preserved, 
is there an optimal choice of initial state, that is, one 
for which the ensuing decoherence is in some sense min- 
imized? To address this question, we must first specify 
what we mean by optimal. Do we wish to minimize aver- 
age decoherence, in which case we would not "see" brief 
but significant drops in the coherence? Alternatively, do 
we wish to minimize the maximum decoherence, in which 
case a significant drop in coherence will make a state ap- 
pear to be a bad choice, even though it may be good 
on average. In the following we will adopt this latter 
criterion. 

It is somewhat easier to use the combinations Hq.i 
rather than H and He. We wish to find the state \T) 
for which the minimum value of \r(t)\ is maximal, in the 
case where several frequencies are present in the sum (J4j) . 
The general case appears difficult to analyze, but one 
might expect that the best choice of |/) is an eigenstate 
of one of the two Hamiltonians and a combination of two 
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eigenstates of the other Hamiltonian. For instance, if 
n* 0) \I) = (II^ + n^ 1} ) \I) = \I), then (0D becomes 

r(t) = e lX ^ (V lA i 4 (I\ \I) + e"^* (I\ H^ |/>) 

and 

|r min | = |(7|n«|7)-(7|n« \I)\. 

This expectation turns out not to be the best choice, in 
general. To see this, we examine the simplest example in 
which no common eigenstate exists, namely, an environ- 
ment consisting of a single qubit with Hq : \ describing its 
interaction with non-parallel magnetic fields of the same 
intensity. Let the precession frequency be oj and let the 
direction of the magnetic fields corresponding to Hq and 
Hi be too = (sin a, 0, cos a) and rhi — (— sin a, 0, cos a), 
with < a < 7r/2, respectively. The initial state can 
be taken to be a spin aligned along a third direction, 
v = (sin 9 cos <f>, sin 9 sin <j>, cos 9) , say. Its evolution ac- 
cording to each of the Hamiltonians is simple: the spin 
precesses around rhi with frequency to so \£i{t)) = \vi(t)) 
where Vi(t) is V rotated about rhi by angle Lot. Then 
r(t) = <er (*)|ei (*)) and \r{t)\ = cos( 7 (i)/2), where j(t) 
is the angle between vo(t) and i>i(t). Thus we would like 
to find the vector v for which the maximum angle be- 
tween i>o(t) and vi(t) as they precess is minimized. If, 
as was conjectured above, we choose v — mo, then it is 
easy to see that the maximum angle is the lesser of 4a 
and 2tt — 4a, reached after half a precession. Although 
it is surprisingly difficult to find this maximum angle for 
an arbitrary v, the optimal choice turns out to depend 
on the angle a as indicated in figure [2j 
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Numerical results 
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Figure 2: Numerical results 

If a < tt/3, it is best to choose v = y. perpendicular 
to too and rhi, resulting in 7 max = 2a (attained after a 
quarter-rotation) and 

| Train | = COS a. 

If a > tt/3, an optimal choice is v collinear with too 
(or mi), resulting in 7 max = 2n — 4a (attained after a 
half-rotation) and 

kminl = C0S(7T - 2a). 



Thus, the optimal value of |r m i n | is close to unity if the 
directions of the two magnetic fields are almost parallel. 
More surprising is the fact that the best choice of initial 
state is far from what one would naively have guessed. 
In the case in which a < tt/3, the optimal choice is v — 
y, which does not correspond to an eigenstate of either 
Hamiltonian. 

To sum up, two regimes emerge. On the one hand, 
for a > tt/3, the system-environment interaction prevails 
whereas the self-evolution of the environment is only a 
perturbation. In this regime, the system and the en- 
vironment play symmetric roles. Reducing decoherence 
in this case boils down to minimizing the entanglement. 
The optimal initial states are thus the pointer states of 
the environment, i.e. the eigenstates of H. On the other 
hand, for a < tt/3, the self-evolution Hamiltonian Hg 
dominates over the interaction Hamiltonian H. However, 
the evolution of the quantum system relies essentially on 
the interaction. Thus, for the evolution of the system, 
H , no matter how small, cannot be considered a pertur- 
bation with respect to Hg. In that case, the evolution of 
the environment is dominated by its own dynamics but 
the impact on the quantum system is mediated by the in- 
teraction Hamiltonian. Hence, finding an analytical cri- 
terion characterizing an initial state of the environment 
that optimally limits the subsequent decoherence of the 
system remains an unresolved challenge. 
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